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We comment on the nature of the ordering transition of a model of equilibrium polydisperse 
rigid rods, on the square lattice, which is reported by Lopez et al. to exhibit random percolation 
criticality in the canonical ensemble, in sharp contrast to (i) our results of Ising criticality for the 
same model in the grand canonical ensemble [Phys. Rev. E 82, 061117 (2010)] and (ii) the absence 
of exponent (s) renormalization for constrained systems with logarithmic specific heat anomalies 
predicted on very general grounds by Fisher [M.E. Fisher, Phys. Rev. 176, 257 (1968)]. 

PACS numbers: 64.60. Cn, 61.20.Gy 



Extensive Grand Canonical Monte Carlo simulations, 
for a model of adsorbed self-assembled rigid rods (S ARR) 
on the square lattice, indicate that the polydisperse rods 
undergo a continuous transition in the two-dimensional 
(2D) Ising class, in line with models of monodisperse rods 
[l|, 0] ■ This finding is in sharp contrast to a previous re- 
sult, based on Canonical Monte Carlo simulations, where 
equilibrium polydispersity was claimed to change the na- 
ture of criticality, from Ising to random percolation . 

In the preceding comment Lopez et al. elaborate on 
this claim to conclude that the criticality of the SARR 
model on the square lattice depends both on polydis- 
persity and on the statistical ensemble. This surprising 
result is based on simulations and normal finite-size scal- 
ing analysis of the SARR model in the canonical and 
grand canonical ensembles. The conclusion was that 
the SARR model exhibits random percolation critical- 
ity (9=1 Potts) when the system is in the canonical 
ensemble, while the criticality is Ising-like {q = 2 Potts) 
when the system is in the grand canonical ensemble. This 
is at odds with very general arguments by Fisher fJ] on 
the absence of exponent renormalization in constrained 
(e.g. fixed density) systems with logarithmic specific heat 
anomalies, as well as with the results of a detailed sim- 
ulation study of Fisher scaling of the 2D Ising magnetic 
lattice-gas, by Ferreira and Prodanescu Fisher has 
also shown that although the universality class of con- 
strained systems, with specific heat anomalies, is that of 
the unconstrained ones there are logarithmic corrections 
to the scaling functions, which may affect the scaling be- 
haviour of reasonably sized systems as shown by Ferreira 
and Prodanescu 

The existence of two universality classes for the SARR 



model, claimed by Lopez et al., is based on the calcula- 
tion of (i) the fourth-order Binder cumulant of the order 
parameter, S, = 1— < > /(3 < >^) at the 
transition, g^, and (ii) the value of the correlation length 
exponent v, obtained by normal scaling data collapse of 
the cumulants, for different system sizes. Both the val- 
ues of g1 and i^, reported by Lopez et al. for the SARR 
model, in the canonical and grand canonical ensembles, 
are different. 

In the canonical simulations of the preceding comment 
Lopez et al. kept the surface coverage constant and var- 
ied the temperature of the system, rather than fixing the 
temperature and varying the coverage f^. As discussed 
below the logarithmic corrections to the normal finite- 
size scaling analysis, arising from the constant density 
constraint apply in both cases. In other words. Fisher 
logarithmic corrections [4, 6] as well as the simpler log- 
arithmic correction suggested by us [l| apply due to the 
constant density constraint (which is also one of the con- 
trol parameters in the canonical ensemble). We stress 
that although Fisher renormalization predicts that the 
critical exponents are unchanged in constrained systems 
with logarithmic specific heat anomalies, as in the SARR 
model on the square lattice, it does predict finite-size 
logarithmic corrections to the scaling functions, which if 
neglected will lead to effective exponents that may differ 
significantly from the true asymptotic exponents of the 
unconstrained system. 

A very careful analysis of the criticality of the Ising 
magnetic lattice gas on the square lattice, in the canoni- 
cal ensemble, was carried out by Ferreira and Prodanescu 
5] and illustrates in detail how the effective exponents 
depend on the scaling analysis of the constrained sys- 
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tern. The authors point out that the values of 174 at the 
intersection of the Binder cumulants for different system 
sizes decrease (slowly) as the system size increases and 
their best estimate for the cumulant at criticality is re- 
ported to be significantly larger than the corresponding 
2D Ising value. Using normal scaling v was found to dif- 
fer from 2D Ising but when Fisher scaling was taken into 
account v was found to approach the 2D Ising value 5] . 
The authors also estimated ^jv and obtained excellent 
agreement with 2D Ising when using Fisher scaling by 
contrast to the value obtained from normal scaling. The 
results for the 2D magnetic Ising gas show clearly that 
when Fisher scaling is taken into account, the effective 
exponents are closer to the values observed in the un- 
constrained system, as expected on theoretical grounds 
0, iSi]. The authors stress that Fisher scaling is not a 
correction to normal scaling but a scaling which devi- 
ates from normal logarithmically, rendering the numeri- 
cal investigation of the criticality of these systems a very 
challenging problem. 

Finite-size-scaling theory asserts that on the critical 
line Tc(/i), g^(V) adopts a non-trivial value, g% indepen- 
dent of the system size L. For a given set of boundary 
conditions, this value of g% is the same for systems in the 
same universality class. In addition, the dependence of 
(74 on the coupling parameter(s), if, in the critical region 
scales asQ: {dg4^/dK) oc L^/'^ . This is what we referred 
to above as normal scaling. Using normal scaling, Lopez 
et al. obtained results for the constrained SARR model, 
consistent with v — 3/4. They also report that the cross- 
ing of gi occurs at g1 ~ 0.638, which is claimed to be the 
value corresponding to the q = 1 Potts universality class 
(random percolation). 

The use of normal scaling for the constrained SARR 
model, leading to Lopez et al. conclusion of percolation 
critical behavior, has to be questioned. Previously [IJ, we 
proposed a simple argument that accounts for the effec- 
tive value of 1/ = 4/3 reported by Lopez et al. Q for the 
constrained SARR model. We indicated that, for large 
systems, there is an additional L/lnL term in the scal- 
ing of the density derivatives compared to field deriva- 
tives. In the range of sizes investigated by [1] L/lnL is 
fitted by: L/lnL ~ aL'^/"' , with i/' ~ 1.291, close to 
1/ = 4/3 of the q = 1 Potts model This logarith- 
mic correction arises from the density constraint and has 
been discussed in much more detail by Fisher and was 
investigated numerically by Ferreira and Prodancscu 5] . 
We note that the simple L/ In L correction [l| is in line 
with Fisher scaling for large systems (see equations (14) 
and (16) of [Ij). Therefore, in what follows we focus on 
the difference between the values of g^ reported by Lopez 
et al., for constrained and unconstrained SARR models. 

Lopez et al., discard the possibility of Ising criticality 
of the constrained SARR model based on the value of 
54, which differs from that of the unconstrained model: 
^ismg ~ 0.611. The results of @ indicate that such an as- 
sumption is far from justified. In normal scaling, appro- 
priate for unconstrained models, finite-size scaling[8j, |9|] 




FIG. 1: Results for the Binder cumulant of the hard square 
lattice model in the grand canonical ensemble for different 
system sizes (See the legends). The filled circle marks the 
critical chemical potential of the HSL model and the corre- 
sponding value of gl for the 2D-Ising universality class (for 
unconstrained systems) with periodic boundary conditions. 



considers the singular part of the appropriate thermo- 
dynamic potential in terms of the thermodynamic fields: 
temperature, T; chemical potential, fi; external fields. 
The corresponding intensive conjugate variables: energy 
per unit volume, density, p; magnetization are the nat- 
ural variables of the constrained models, where Fisher 
scaling applies [3, Hj. The SARR model on the square 
lattice may be described as a symmetric binary mixture, 
where a species corresponds to a given orientation. The 
relevant thermodynamic fields are then T and /z. Within 
this Grand Canonical description of the SARR model 
(completed by taking the volume as the extensive vari- 
able that defines the system size) , normal scaling theory 
applies. Of course, one can invesigate the criticality of 
the model in other ensembles but then the appropriate 
scaling theory must be used [1, [13] ■ 

In order to check the effect of the density constraint 
on g^ in systems where {dp/dii)T diverges at the critical 
point we consider the behavior of the hard square lattice 
(HSL) model [ll|. The HSL is an athermal model (an 
occupied site excludes occupation of its nearest neighbor 
sites) defined on the square lattice and exhibits a contin- 
uous order-disorder transition: at high densities particles 
occupy preferentially one of the two sublattices. The or- 
der parameter is defined as: 5 — |iVi — iV2|/L^, where 
Ni is the number of occupied sites in sublattice i, and 
L^ is the number of lattice sites. The transition of the 
HSL model is in the 2D Ising class and both the chemical 
potential and the density at the critical point are known 
with high accuracy [11]. 

We simulated the transition of the HSL model using a 
multicanonical sampling procedure [ll, [T2| - [l^ that allows 
results in the canonical and grand canonical ensembles to 
be obtained simultaneously. In figures [T] and [2] we illus- 
trate the results for 174 (/3^) and (?4(p), in the grand canon- 
ical (unconstrained) and canonical (constrained) ensem- 
bles respectively. We find that the results in the grand 
canonical ensemble are fully consistent with the expected 
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FIG. 2: Results for the Binder cumulant of the hard square 
lattice model in the canonical ensemble. The filled circle 
marks the critical density of the HSL model and the corre- 
sponding value of for the 2D-Ising universality class (un- 
constrained systems) with periodic boundary conditions, po 
is the density at maximum lattice occupancy. 

2D Ising behavior. The curves g4{/3fi) for different sys- 
tem sizes cross (within error bars) at the the expected 
value {Pfj,CT g^^"'^). However, in the canonical ensemble, 
the crossings occurs at a density slightly larger than p^, 
(this could be a finite-size effect), while the crossing of 34 
decreases slowly as the lattice size increases, in line with 
the results reported for the Ising lattice gas model 
More importantly, the results suggest that the universal 
value of (74 for the constrained system may differ from 



the 2D Ising value for the unconstrained system, g^^^^ . 
Incidentally, the crossings of 174 in the canonical ensem- 
ble, occur at values close to the value reported by Lopez 
et al. as the universal value of the cumulant for the q = 1 
Potts criticality. 

We conclude that the dependence of the universality 
class of the SARR model on the statistical ensemble, re- 
ported by Lopez et al., is very likely the result of inad- 
equate use of normal scaling to investigate the critical 
properties of the constrained (constant density) system. 
A full analysis following the lead of Ferreira and Pro- 
danescu @ seems to be called for but it is clearly outside 
the scope of this Reply. 
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